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INTRODUCTION 

The  fiber/matrix  interface  has  a  significant  influence  on  the  structural  integrity  of  a  fibrous 
composite.  For  example,  the  stronger  the  interface  bond,  the  higher  the  composite's  static  strength. 
However,  a  strong  interface  bond  also  yields  a  brittle  and  notch-sensitive  composite*  and,  on  the  other 
hand,  a  weaker  interface  results  in  a  higher  fracture  toughness  composite.2,3 

Conventional  methods  for  predicting  the  stress  at  the  fiber/matrix  interface  can  be  divided  into  three 
categories.  These  are:  shear  lag  theory,4,5,4  two-dimensional,  finite-element  analysis, 7-4  and  three- 
dimensional,  finite-element  analysis.6  These  methods  all  assume  zero  interface  thickness  and  uniform, 
homogeneous  matrix  properties. 

Recently,  several  researchers  have  suggested  that  the  volume  of  material  immediately  surrounding 
the  fiber  is  significantly  different  from  the  bulk  matrix.10,11, 12,13  This  volume  of  material  is  commonly 
referred  to  as  the  interphase.  Drzal10  has  suggested  that  this  material  may  be  more  rigid  than  the  bulk 
matrix.  Piggott11  has  found  that  to  explain  the  Young's  moduli  of  short-fiber  composites,  the  interphase 
must  have  a  very  low  modulus,  i.e.,  much  softer  than  the  bulk  matrix. 

No  direct  evidence  has  been  found  for  the  presence  of  an  interphase  in  organic  matrix  composites 
nor  are  its  properties  and  dimensions  known.  The  objective  of  this  study  is  to  prove  the  existence  of 
interphase  and  to  determine  interphase  elastic  properties,  if  possible.  Included  in  this  work  is  an 
investigation  to  determine  the  influence  of  interphases  on  stress  at  the  fiber /matrix  interfaces  on 
displacements  and  on  fracture  toughnesses  of  fibrous  composites. 


SHEAR  LAG  THEORY  WITH  DISTINCT  INTERPHASE 

To  include  the  interphase  into  a  composite-material  structural  analysis  model,  material  properties 
and  the  thickness  of  the  interphase  must  be  known  in  advance.  Experimental  data  on  interphase 
properties  does  not  exist.  One  way  to  deduce  the  interphase  characteristics  is  to  assume  initial  values 
for  the  interphase  and  then  iterate  the  properties  and  thickness  until  analytical  results  converge  to  match 
the  corresponding  experimentally  observable  results. 

Since  the  thickness  of  interphase  can  be  very  small  compared  with  the  fiber  diameter,  a  very  large 
number  of  finite  elements  would  be  needed  if  the  finite  element  technique  was  used  for  the  iteration 
procedure.  This  would  be  time-consuming  and  costly.  A  closed-form  solution,  based  upon  shear-lag 
theory,  was  employed  instead. 

In  the  following,  a  shear  lag  analysis,  which  includes  and  interphase  region  between  the  fiber  and 
bulk  matrix,  has  been  developed.  Mandell’s  microdebonding  test  data1  was  used  to  determine  the 
interphase  material  properties  and  thickness. 

THEORY 

Mandell  s  microdebonding  test,1  Figure  1,  was  modeled  as  shown  in  Figure  2.  Four  types  of 
materials  were  included  in  this  model;  namely,  fiber,  interphase,  matrix,  and  composite.  Note  that  a  load 
is  applied  at  the  end  of  a  fiber,  and  that  the  material  outside  the  matrix  was  modeled  as  equivalent 
smeared  composite  material. 

By  employing  the  shear  lag  theory,  the  following  assumptions  were  made: 

(1 )  Axial  load  is  carried  by  the  fiber  alone,  while  the  interphase  and  matrix  carry  the  shear  load  only. 
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LOAD 


Figure  1.  Schematic  of  Loading  for  Microdebonding  Test. 
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Figure  2.  Modeling  of  Fiber/Matrix  Interlace 
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(2)  Fixed  boundary  conditions  apply  at  the  boundary  between  the  matrix  and  equivalent  composite. 

The  freebody  diagrams  of  the  loaded  fiber  and  adjacent  interphase  and  matrix  are  shown  in  Figures 
3a  and  3b.  The  equilibrium  equations  are  as  follows: 


■~+2nrtxi  =  0 
dz  '  ' 


(i) 


~+2nrx  ~  0 
dz  m  m 


(2) 


From  equations  (1)  and  (2),  we  have: 

=  rmXm  =  rx  (3) 

where 

F  =  fiber  force  at  z 

rf,  r,  =  radius  of  fiber  and  interphase,  respectively 

r,,  r„  =  shear  stress  of  interphase  and  matrix,  respectively. 

From  the  theory  of  elasticity  and  Figure  4: 


fF<  r<  rt 
dW_  _  _x_ 

dr  =  Gi  (4) 

where  r  is  the  shear  stress  in  the  interphase,  and  G,  is  the  shear  modulus  of  interphase. 
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From  equation  (3): 


substituting  equation  (5)  into  (4)  and  integrating,  we  have: 


w  =  —pr~ln(r)+  C 

Gi 


r  =  ry  w  =  wj 


r  =  rlt  vt>  =  v>2 


Applying  these  boundary  conditions  to  equation  (6),  we  have: 


(w2- wjjG, 


r,/n(— ) 

7 


By  the  same  procedure  and  noting  that  at  r=rm,  w  =  0  we  have: 


-w2Gn 


-  d  w. 

F  =  Efnrf2 - 

7  7  dz 
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substituting  equations  (7)  through  (9)  into  equation  (1)  and  (2)  we  have: 


d2wx 
dz 2 


2  Gi 


1  1  7 


■(h'1_w2)  =  0 


d2wx 
dz 2 


2G. 


rfEfln(-j-) 


■w2 


=  0 


After  some  algebraic  manipulations  of  equations  (10)  and  (11),  we  have: 


where 


d2wx 
dz 2 


-a2wl  =  0 


2Gj 

Err2ln(-~~) 

rf 

Giln  (——) 

1+ - T~ 

7 


The  solution  of  equation  (12)  is. 

w  i=  C  \cosh  ( az)+C  2sinh  ( a  z ) 


(10) 


(ID 


(12) 


(13) 
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z  =  0,  F  =  -Pq  =  -  7t  rffS  -  Ejitrf 


2  =  1,  Wj  =  0 


where  5  is  the  average  applied  stress  at  the  fiber  end. 

Applying  the  boundary  conditions  (l.e.,  equation  (14))  to  equation  (13),  we  have: 


wj  =  -  — — (-tanh(al)cosh  (az)+ sinh(az)) 


w2  is  obtained  from  equations  (11)  and  (15)  such  that: 


r?EFa2ln(-?-) 


w2  = 


From  equations  (7),  (15),  and  (16),  the  interphase  shear  stress  at  r=r,  is  written  as: 

T  =  ^  Gm  rf  (-tanh(al)cosh(az)+sinh(az)) 

*  V  2 Et  r.-  rm  ~~ 

\  V,-  r/  r, 

From  equation  (3),  the  interphase  shear  stress  at  r=ry  is  written  as: 


T*  =  VI 


(- tank (al) cosh (a.z)+sinh(az} 

V  — ' t 


The  maximum  interphase  shear  stress  occurs  at: 


_ron/i(a/) 
Gm  Ti 


In  (-—■)+ In 


Gi  'r. 


Note  that  equation  (19)  contains  two  unknowns  to  be  determined  (namely,  G,  and  r). 
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CORRELATION  WITH  MANDELL’S  MICRODEBONDING  TESTS 

DETERMINATION  OF  INTERPHASE  PROPERTIES 

We  now  postulate  that  debonding  of  fiber  from  the  matrix  Is  due  to  interphase  shear  failure.  Based 
on  this  assumption,  the  debonding  criterion  for  the  fiber/matrix  interface  is  written  as  follows: 


where 

Tq  =  interphase  shear  strength 

s0  =  average  applied  stress  at  the  fiber  end  which  causes 

fiber-matrix  debonding  (referreo  to  as  the  debonding  stress). 

Let 

r-  =  rf  +  t. 
rm  =  rf+tm 

where  t,  is  interphase  thickness  and  tm  is  the  thickness  combining  interphase  and  matrix. 
Equation  (20)  becomes: 


From  equation  (21 )  we  conclude  that: 


By  selecting  tm/d, =0.4  and  tjd, = 1 .0  to  correlate  the  test  data,1  the  following  results  were  obtained 
(for  detailed  calculations,  refer  to  Appendix  A). 
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For  S-glass/epoxy, 


G  2  r 

(•7T--1)/«(1+-Ti)=  1.8446 

Gi  df 

x0  =  39  MPa 


For  graphite/epoxy 


(23) 


G  2 1 

(-~-l)/«(l+-~)  =  2.0672 
G,  A/ 


T0  =  27  AfPa 


(24) 


Note  that  Gi  and  t,  cannot  be  determined  separately  from  equations  (23  and  (24),  but  if  f,  can  be 
determined  from  an  experiment,  then  G,  can  be  determined.  Tables  1  and  2  list  the  interphase 
properties  of  S-glass/epoxy  and  graphite/epoxy,  respectively,  for  various  interphase  thickness  values 
and  show  the  interphase  shear  modulus  to  be  softer  than  bulk  matrix  for  both  S-glass/epoxy  and 
graphite/epoxy  composite  materials. 

CORRELATION  WITH  MICRODEBONDING  TEST  RESULTS 

After  the  interphase  shear  strength  has  been  determined,  the  debonding  stress  d0  is  determined  as 
follows: 


O0  =  VTt0 


4J'! 


2  h 


It, 


(~7T —  l)ln(l+~~r-)+ln(l+-r  ) 
Gi  df  df 


(25) 


Making  use  of  equations  (23)  and  (24),  and  noted  that  we  have  for  S-glass/epoxy: 


=  T0 r\j  2 1 .8446+  In  ( 1+ 


(26) 


and  for  graphite/epoxy: 


=  t0 yj  2-^-^j  (2.0672+  /«d+-^L) 


(27) 


Equations  (26)  and  (27)  are  plotted  in  Figures  5  and  6.  The  correlations  with  the  microdebonding 
test  results  are  excellent. 
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Table  1.  Interphase  Properties  of  S-glass  Epoxy 


Interphase  Shear  Strength 

=  39  MPA 

tj/df 

t|  (nm) 

Gm''Gi 

.001 

10 

922 

.005 

50 

184 

.01 

100 

92 

.05 

n 

18 

.1 

1000 

9.0 

Table  2.  Interphase  Properties  of  Graphite/Epoxy 


Interphase  Shear  Strength  =  27  MPA 


tj/df 

tj  (nm) 

.0014 

_ 

10 

724 

.0071 

50 

145 

.0143 

100 

74 

.0714 

m 

15 

.143 

1000 

7 

0.9 
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Figure  5.  Variation  of  Debond  Stress  vs.  t^/d^  for  S-glass/Epoxy. 
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Figure  6.  Variation  of  Debond  Stress  vs.  for  Graphite/Epoxy. 
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INFLUENCE  OF  INTERPHASE  PROPERTIES  ON  INTERPHASE 
SHEAR  AND  FRACTURE  TOUGHNESS  OF  COMPOSITES 

INTERPHASE  SHEAR  STRESS 

Equation  (19)  can  be  rewritten  as: 


Equation  (30)  is  plotted  in  Figures  7  and  8  for  Ks  vs.  Gm/G,  for  various  ratios.  We  can  conclude  that: 

(1)  The  lower  the  Gm/G,  ratio,  the  higher  the  stress  concentration  and  the  more  likely  debonding. 

(2)  The  thinner  the  interphase  thickness,  tit  the  higher  the  stress  concentration,  again  making 
debonding  easier. 

MODE  I  FRACTURE  TOUGHNESS 

Murphy,3  etc.,  showed  that  the  Mode  I  fracture  energy  of  unidirectional  composites  can  be 
expressed  as  follows: 

(G/)c  *  0-Vf)GM  +  Vf(GF+GINF)  (31) 

where 

(Gi)c  = 


Mode  i  fracture  energy  of  composite 


Figure  7.  Effect  of  Interphase  Shear  Modulus  on  Maximum  Stress  for  tm/d^,=0.1. 
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Figure  8.  Effect  of  Interphase  Shear  Modulus  on  Maximum  Interfacial  Shear  Stress  for  t _/d,  =  1.0. 
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Gm  =  Mode  I  fracture  energy  of  matrix 

Gf  =  Mode  I  fracture  energy  of  fiber 

Vf  =  volume  fraction  of  fiber 

G,nf  =  Mode  I  fracture  energy  associated  with  the  decohesion  of  the  fiber  matrix  interface. 

Murphy  further  showed  that: 


dfah  ,  ,  1  lD  2  Et 

0,„F= 

4V£/  2  l*  °/« 


(32) 


where 

°U 

I a 


ultimate  stress  of  fiber 

(d//4r/)  (af{J  -  50) 

debonding  stress 

fiber  pull-out  length 

post  debonding  friction  shear  stress 

debonding  fracture  ratio 


(33) 


Examining  equations  (31)  and  (32),  it  can  be  seen  that  the  smaller  the  g  (i.e.,  weaker  interface),  the 
larger  the  GINC  and  hence,  the  tougher  the  composite.  On  the  other  hand,  the  larger,  the  g  (i.e., 
stronger  interlace),  the  smaller  the  G,Nr  yielding  a  more  brittle  composite.  If  g=  1 ,  (i.e.,  the  fiber  breaks 
with  no  debonding  between  fiber  and  matrix),  equation  (31)  reduces  to: 

<Vc  ■  <1-V,)  CM  -  vfCF 


Murphy3  showed  that  for  high-volume  fractions  the  Mode  I  fracture  energy,  (G,)c,  of  composites  can 
range  over  five  orders  of  magnitude,  while  the  debonding  fracture  ratio,  g,  ranges  from  zero  to  one. 

From  equations  (29)  and  (33),  g  can  be  expressed  as  follows: 


Ofu  Ks 


(34) 


18 


N  ADC-90001 -60 


To  demonstrate  how  interphase  properties  will  effect  the  Mode  I  fracture  toughness  of  the 
composite,  equation  (32)  is  rewritten  as  follows: 


ginf 

ClNF 


1  -g+g2-g3+lr(^-)2(l-g)2 
z  ld 


(35) 


where 


dfGfu 
ClNF  =  4X/£/ 


For  simplicity,  assume  /p=0  (i.e.,  no  fiber  pull-out  occurs  during  fracture  of  composite).  Equation 
(35)  then  becomes: 


G1NF 

CINF 


l-g+g2~g 3 


(36) 


Assuming  that  the  interphase  strength,  t0,  is  constant  for  various  interphase  thicknesses  and  making 
use  of  equation  (34),  equation  (36)  is  plotted  as  shown  in  Figures  9  and  10.  The  detailed  calculation 
is  shown  in  Appendix  B. 

From  the  plots,  it  can  be  seen  that  when  the  interphase  property,  G„  is  closer  to  the  bulk  matrix, 
Gm,  the  composite  material  will  have  higher  fracture  toughness.  Note  that  the  thinner  the  interphase 
thickness,  f„  the  tougher  the  composite  material. 


FINITE  ELEMENT  SIMULATION 


MICRODEBONDING  PROBLEM 

An  axisymmetric  finite  element  model  that  includes  the  interphase  was  developed.  Interphase 
properties  were  obtained  by  either  using  the  values  contained  in  Tables  1  and  2  or  those  derived  from 
using  equation  (23)  for  S-glass/epoxy  and  equation  (24)  for  graphite/epoxy  for  various  interphase 
thicknesses.  Only  those  results  for  graphite/epoxy  are  presented.  Figures  1 1  and  1 2  are  typical  models 
used  to  analyze  the  microdebonding  problem.  The  NASTRAN  computer  code  was  used  to  perform  the 
analysis.  Note  that  the  top  surface  and  side  surface  corresponding  to  outside  radius  were  assumed 
fixed.  A  1  N  force  was  applied  at  the  center  of  the  cylindrical  composite  disc. 

Figure  13  shows  the  typical  interphase  shear  distribution  along  the  axis.  The  maximum  shear  stress 
occurs  at  approximately  0.4  fiber  diameter  below  the  surface.  Shear  lag  theory  positions  the  maximum 
shear  at  the  surface.  Note  that  the  magnitude  of  the  interface  shear  for  composite  without  an  interphase 
(i.e.,  Gm/G,  =  l.O)  is  much  higher  than  for  a  composite  with  an  interphase.  This  result  agrees  with  the 
conclusion  from  the  shear  lag  analysis  in  the  foregoing  section. 
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Figure  9.  Effect  of  Interphase  on  Opening  Mode  Fracture  Energy  of  a  Unidirectional  Fiber  Composite  for  tm/d^=0.1. 
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Finite  Element  Model  for  Axisymmetric  Solid. 
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Figure  14  shows  the  prediction  of  debonding  stress  d0,  for  tm/df-0A  and  Xm/df=t.O,  where  data 
for  tjdf= 0.4  was  used  to  determine  the  interphase  shear  strength.  We  see  that  the  model,  which 
includes  the  interphase,  provides  a  better  solution  than  those  without  the  interphase. 

Also,  note  that  interphase  shear  strength  predicted  by  finite  element  analysis,  23.9  mpa,  is  close  to 
that  predicted  by  the  shear  lag  theory.  (For  detailed  calculation,  refer  to  Appendix  A.  This  interphase 
shear  strength  value  is  different  from  the  values  of  interfacial  shear  strength  obtained  by  other 
investigators.' 

Figure  15  shows  the  fiber  end  displacement  in  the  axial  direction  under  a  1.0  N  axial  load  applied 
at  the  fiber  end.  The  correlation  between  the  finite  element  analysis  and  new  shear  lag  theory  is  very 
good,  validating  the  adequacy  of  the  shear  lag  theory. 


INTERPHASE  EFFECT  ON  OUT-OF-PLANE  DISPLACEMENT 

Finite  element  calculations  were  performed  to  determine  the  effect  of  the  interphase  on  the  out-of¬ 
plane  displacement  using  the  model  shown  in  Figure  16.  This  model  corresponds  to  the  experimental 
geometry  used  for  tensioned  fiber  tests  in  reference  14.  Figures  17  and  18  show  the  predicted  out-of¬ 
plane  displacements  with  and  without  interphase  for  Gm/G=25  and  Gm/G,=74  respectively.  It  can  be 
seen  that  a  significant  displacement  gradient  occurs  in  the  vicinity  of  interphase  region.  Thus,  the 
interphase  region  may  be  defined  as  a  region  with  a  very  high  displacement  gradient  allowing  the 
dimension  of  interphase  to  be  determined  through  direct  measurement  of  the  length  of  this  high 
displacement  gradient  region. 

The  results  shown  in  Figures  17  and  18  were  adjusted  for  a  0.4%  fiber  strain  load  and  correlated 
with  test  results'5  as  shown  in  Figure  19.  Interphase  properties  with  GJG,- 25  and  t^O.Um  match  the 
test  data.  However,  if  this  method  is  not  believed  adequate  to  determine  the  exact  interphase  thickness, 
then  the  radial  displacement  field  must  also  be  measured  in  addition  to  the  axial  displacement  field  so 
that  the  test  results  can  be  correlated  to  determine  separate  values  for  the  interphase  thickness  and 
modulus.  Note  that  the  softer  the  interphase,  the  greater  the  displacement  gradient  in  the  interphase 
region. 


DISCUSSION  AND  CONCLUSION 

The  ratio  between  interphase  shear  modulus  and  thickness  can  be  established  from  shear  lag 
analysis,  but  not  specific  values  for  G(  or  T ,  separately,  if  the  interphase  thickness  can  be  determined 
from  an  experimental  measurement,  then  equation  (23)  or  equation  (24)  can  be  used  to  calculate  the 
interphase  shear  modulus  for  S-glass/epoxy  and  graphite/epoxy,  respectively.  If  the  interphase 
thickness  is  not  available  experimentally,  then  an  iteration  procedure  must  be  performed  using  equations 
(23)  and  (24)  as  starting  values  for  interphase  shear  modulus  and  thickness,  iterate  G,  and  T,  until  the 
analytical  results  (i.e.,  out-of-plane  displacement  and  radial  displacement)  converge  to  the  test  results. 
More  studies  need  to  be  pursued  of  this  approach  to  compute  G;  or  Tr 

In  this  report  it  is  assumed  that  the  interphase  is  an  axisymmetric  region  of  uniform  modulus  and 
fixed  radial  dimension.  We  recognize  that  the  true  interphase  region  has  radially  varying  properties  and 
can  also  exhibit  circumferential  property  gradients. 

For  a  tougher  composite,  equation  (32)  and  Figures  9  and  10  suggest  that  the  composite  material 
must  possess: 
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1 .  Larger  fiber  diameter  and  fiber  ultimate  stress 

2.  Smaller  post  debond  frictional  shear  stress  and  fiber  Young’s  modulus 

3.  Smaller  interphase  shear  strength 

4.  Smaller  interphase  thickness  and  Gm/G,  ratio. 

Although  interphase  shear  strength,  r0,  and  post  debond  frictional  shear  stress  vary  with  interphase 
thickness.  Figures  9  and  10  give  approximate  guidelines  for  obtaining  tougher  composite  materials. 
Note  that  the  weaker  interface  (i.e.,  small  interphase  shear  strength)  provides  tougher  fracture 
resistance,  but  may  affect  the  strength  of  composite,  especially  short  fiber  composites. 

The  following  conclusions  are  drawn  from  this  study: 

1 .  By  including  the  interphase,  both  shear  lag  analysis  and  finite  element  analysis  provide  much 
better  predictions  of  material  response. 

2.  The  interphase  exists  and  is  softer  than  the  matrix  for  the  composites  used  by  Mandell  in  his 
microdebonding  tests  (i.e.,  S-glass/epoxy  and  graphite/epoxy). 

3.  The  properties  and  thickness  of  the  interphase  have  significant  influence  on  the  interface  stress, 
displacement,  and  fracture  toughness  of  fibrous  composites. 
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With  and  Without  Interphase  for  Graphite/Epoxy. 
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:ial  Displacement  of  Fiber  End  Under  1  N  Load  Applied  at  Fiber  End. 
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Figure  18.  Interphase  Effect  on  Displacement  at  Top  Surface 
(P  =  -0.017N,  tj  =0.1  <«m),  G^G,-  =74). 


VERTICA 


RADIAL  POSITION  ( fim ) 


Figure  19.  Correlation  of  Finite  Element  Analysis  and  Out-of-Plane  Displacement 
Test  Data  for  0.4%  Fiber  Strain  Load  With  t  =  0.ln  m. 


N  ADC-90001 -60 


REFERENCES 

1.  Mandell,  J.F.,  et  al.  "Modified  Micro-debonding  Test  for  Direct  In-Situ  Fiber/Matrix  Bond  Strength 
Determination  in  Fiber  Composites,"  Composite  Materials  (Seventh  Conference).  ASTM  STP  893. 

2.  Chamis,  C.C.  "Mechanics  of  Load  Transfer  at  the  Fiber/Matrix  Interface,"  NASA  TN  D-6588,  February 
1972. 

3.  Murphy,  M.C.  and  Outwater,  J.O.  "The  Toughness  of  Reinforced  Plastics,"  28th  Annual  Technical 
Conference  1973.  Reinforced  Plastic/Composites  Institute.  The  Society  of  Plastics  Industry,  Inc. 

4.  Greszczuk,  L.B.  Theoretical  Studies  of  the  Mechanics  of  the  Fiber-Matrix  Interface  in  Composites," 
ASTM  STP  452,  American  Society  for  Testing  and  Materials,  Philadelphia,  p.  42,  1969. 

5.  Cox,  H.L.  Br.  J.  Appl.  Phys.  16,  p.  199  (1952). 

6.  Piggott,  M.R.  “Load  Bearing  Fibre  Composites"  Pergamon,  Oxford  (1980). 

7.  Remedios,  N.C.  and  Wood,  W.G.  "Stress  Transfer  from  a  Loaded  Matrix  to  a  Single  Fiber,"  J. 
Composite  Material,  Vol.  2,  1968,  pp.  517-520. 

8  Carraron,  A.S.  and  McGarry,  F.J.  "Matrix  and  Interface  Stresses  in  a  Discontinuous  Fiber  Composite 
Model,"  J.  Composite  Mat.  Vol  2,  No.  2,  April  1968,  p.  222. 

9.  Harner,  Juan,  Ashbaugh,  Noel,  et  al.  "Investigation  of  Micromechanical  Behavior  of  Fiber  Reinforced 
Plastics,"  Wittaker  Corp.  (USAAVLABS-TR-67-66,  AD-667901),  February  1968. 

10  Drzal,  Lawrence  T.  Tough  Composite  Materials:  Recent  Developments,"  Published  by  Noyes  Pub., 
Park  Ridge.  NJ,  1985,  pp.  207-222. 

11.  Piggott,  M  R.  Polymer  Composites  8  (5),  1987,  pp  291-296. 

12.  Garton,  Andrew;  Daly;  John  H.  Polymer  Composites  6  (4),  1985. 

13.  Shelton,  C.G.  and  Marks,  P.R.  Journal  of  Material  Science  Letters  7  (6)  1988,  pp.  673-675. 

14.  Williams,  J.G.,  Donnellan,  M.E.,  et  al.  "Properties  of  the  Interphase  in  Organic  Matrix  Composites," 
Journal  of  Material  Science,  January  1990. 

15.  James,  M.R.,  Morris,  W.L.,  and  Dadkhah,  M.  "Determination  of  interfacial  Properties,"  Interim  report 
for  period  1  March  through  15  June  1989.  Contract  No.  N00014-89-C-0083,  Rockwell  International 
Science  Center,  June  1989. 


33 


N  ADC-90001 -60 


APPENDIX  A:  DETERMINATION  OF  INTERPHASE  PROPERTIES 
FOR  S-GLASS/EPOXY  AND  GRAPHITE/EPOXY 

NEW  SHEAR  LAG  THEORY 

INTERPHASE  PROPERTIES 


From  equation  (22),  we  have: 


q  tanh  (a  t) 


1/2 


Const 


(A-1) 


Substituting  tjdf=  1.0  and  tm/df=0A  into  equation  (A-1),  and  note  that  tanh  (a  /)= i.o,  we  have: 

°i.o  °o  A 


(A  +  1 ,0986)1/2  (A  ♦  0.5878)1/2 


(A-2) 


where 


°o.*’  5,.o 


debonding  stress  of  fiber/matrix  Interface  corresponding  to  tjdf- 1.0  and 
tJdf=QA,  respectively 


A 


(A-3) 


Equation  (A-2)  can  be  rearranged  as  follows: 


(Gm  W 
—^-1  In 

1 

(1.0986  -  0.5878  B) 

U  J  l 

'J 

(B  -  1) 

(A-4) 


where 


B  =  (3,Joot)! 

For  S-glass/epoxy,  from  reference  1 : 


(A-5) 


°o.« 


B  -  1.21 


A-1 
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equation  (A-4)  becomes: 


For  graphite/epoxy,  from  reference  1 : 

■!=  -  1.092 

°0.4 

B- 1.1 924 

equation  (A-4)  becomes: 

-A) -2.0672 


INTERPHASE-SHEAR  STRENGTH 


From  equation  (28),  the  maximum  interfacial  shear  stress  can  be  written  as  follows: 

1-1/2 


y/2  \  5 


(A-6) 


For  S-glass  fiber  Gm=1.07  GPA,  E/=86  GPA  from  reference  1,  for  tm/df =0.4,  5=777  MPA, 
equation  (A-6)  becomes: 

r0  =  39  MPA 

For  G,/Gp,  Gm  =  1.21  GPA,  E/=241  GPA,  from  reference  1,  for  tm/df =0.4,  5=883  MPA,  equation 
(A-6)  becomes: 

r0  =  27  MPA 

OUT-OF-DISPLACEMENT  OF  FIBER 


From  shear  lag  theory,  the  fiber  axial  displacement  at  z=0  can  be  written  as: 

o 


w  - 


E,a 


(A-7) 


A-2 
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where 


mJ& 
r*  N 


l  °i  )  UJ  UJ, 


-1/2 


For  graphite/epoxy  an  P=1  N,  equation  (A-7)  becomes: 

a  25.984  1  0.1078 


w  - 


Et  s  241  a 


a  -  0.0286 


2.0672 


■(■  •  *) 


■1/2 


DEBONDING  STRESSES 
Graphite/Epoxy 

From  equation  (A-1),  we  have: 


'1.0 


1 2.0672  + 

After  simplification,  we  have: 


2f 

In  1  +  — -p 


-1/2  {  2.0672  +  1.0986}1/z 


o  -  542.36 


2.0672  +  In 


R) 


1/2 


o  vs.  tm/d/  is  tabulated  as  shown  in  Table  A-1. 
S-glass/Epoxy 

From  equation  (A-1)  we  have: 


'1.0 


2f 

In  1  ♦ 

df) 


|1.8446  + 

After  simplification  we  have: 

o  -  498.375 


1/2  { 1.8446  +  1.0986}1* 


1.8446 


1/2 


A-3 


(A-8) 

(A-9) 

(A-10) 
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tjdf  is  tabulated  as  shown  in  Table  A-2. 


Table  A-1 .  5  vs.  t^/d/  for  Graphite/Epoxy 


U/d, 

3^  (MPA) 

a  fMPAI 

Ep/o 

0.1 

773.6 

813.45 

5.2 

0.2 

826.8 

840.86 

1.7 

0.3 

859.7 

863.9 

0.49 

0.4 

883.0 

883.0 

0. 

0.5 

902.9 

901.09 

-0.2 

0.6 

918.8 

916.51 

-0.25 

0.7 

932.5 

930.37 

-0.23 

0.8 

944.5 

942.94 

-0.17 

0.9 

955.3 

954.4 

-0.09 

1.0 

965.0 

965.0 

0. 

Ud, 

Table  A-2.  a  vs. 

3™  (MPA) 

L/d  /  for  S-gtass/Epoxy. 

5  /MPAI 

Eo^ 

0.1 

673.0 

709.5 

5.4 

0.2 

723.0 

736.0 

1.8 

0.3 

754.0 

758.2 

.56 

0.4 

777.0 

777.0 

0. 

0.5 

795.5 

793.92 

-0.2 

0.6 

810.8 

808.7 

-0.26 

0.7 

823.9 

821.95 

-0.24 

0.8 

835.4 

833.95 

-0.17 

0.9 

845.7 

844.92 

-0.09 

1.0 

855.0 

855.0 

0. 

A-4 
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FINITE  ELEMENT  ANALYSIS 

From  finite  element  analyst,  for  center  fiber  subjected  to  5=25,984.0  MPA  (i.e.,  P=1  N),  the 
computer  outputs  for  interfacial  shear,  r,  and  out-of-piane  displacement,  w,  are  shown  in  Table  A-3. 


Table  A-3. 

Finite  Element  Results  for  P= 

1  N  (t,=.1|im). 

Ln/d, 

w 

T 

(anil 

(MPA) 

0.10 

5.62 

775.0 

0.35 

6.23 

710.7 

1.0 

6.86 

641.0 

INTERPHASE  SHEAR  STRENGTH 


Define  the  interphase  shear  stress  at  5=  883  MPA,  which  corresponds  to  the  debonding  stress  of 
composite  with  tjdf =0.4  as  interphase  shear  strength,  then  r0= 704/25984*883  =23.9  MPA. 

DEBONDING  STRESS  5  vs.  tjdf 


The  debonding  stress  of  composite  for  various  tjdf  can  be  written  as  follows: 


o 


2§§§!x23.9 


(A-11) 


where 

r*  =  computer  output  of  interphase  shear  stress  corresponding  to  5=25984  MPA. 

Making  use  of  data  from  Table  A-3,  we  have  the  comparison  of  finite  element  results  with  the  test 
results  as  shown  in  Table  A-4. 


Table  A-4.  Comparison  of  Debonding  Stress 
Between  Test  and  Finite  Element  Analysis. 


T* 

5 

Ln/d, 

(MPA) 

(MPA) 

MPA 

0.1 

773.6 

775.0 

801.3 

0.35 

874.2 

710.7 

874.2 

1.0 

965.0 

641.0 

969.2 
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OUT-OF-PLANE  DISPLACEMENT 

The  comparison  of  oi* -"f -phase  displacement  for  shear  lag  theory  and  finite  element  analysis  is 
shown  in  Table  A-5.  Shear  lag  solution  Is  based  on  equations  (A-9)  and  (A-10). 


Table  A-5.  Comparison  of  Out-of-Plane  Displacement  Between 
Shear  Lag  Theory  and  Finite  Element  Analysis. 


w(pm) 

SHEAR 

wOtm) 

Ln/d, 

LAG 

E£ A 

0.1 

5.653 

5.62 

0.2 

5.84 

- 

0.3 

6.004 

- 

0.4 

6.075 

6.23 

0.5 

6.2623 

- 

0.6 

6.37 

- 

0.7 

6.47 

- 

0.8 

6.55 

- 

0.9 

6.63 

- 

1.0 

6.71 

6.86 

A-6 


NADC-90001-60 


APPENDIX  B:  CALCULATION  OF  INTERPHASE  SHEAR  STRESS 
CONCENTRATION  FACTOR  AND  MODE  I  FRACTURE 
ENERGY  OF  UNIDIRECTIONAL  COMPOSITES 

INTERPHASE  SHEAR  STRESS  CONCENTRATION  FACTOR,  Kg 

From  equation  (30)  we  have: 


Equation  (B-1)  is  used  to  calculate  Kt  for  fm/d,=0.1  and  fm/cf,= 1.0,  respectively,  as  tabulated  in 
Tables  B-1  and  B-2. 


Table  B-1.  K,  vs.  Gm/G,  for  Various  t,/d,  (t„,/d(=0.l). 


K. 


Gm/G, 

t,/d,=0.001 

t,/d,= 0.005 

t,/df  =  0.01 

1 

1.66 

1.66 

1.66 

10 

1.58 

1.36 

1.17 

102 

1.15 

0.65 

0.48 

103 

0.48 

0.22 

0.16 

104 

0.16 

0.1 

0.05 

105 

0.05 

0.02 

0.02 

Table  B-2.  K,  vs.  Gm/G,  for  Various  tyd,  (tm/d,=i  .0). 

K, 


GrrVG, 

ti/d, =0.001 

t,/d,= 0.005 

t,/d,=0.01 

1 

0.675 

0.675 

0.675 

10 

0.669 

0.649 

0.626 

10* 

0.621 

0.49 

0.404 

103 

0.402 

0.21 

0.155 

104 

0.15 

0.07 

0.05 

105 

0.05 

0.02 

0.02 
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MODE  1  FRACTURE  ENERGY  OF  UNIDIRECTIONAL  COMPOSITE 
From  equation  (29): 


The  debonding  stress  of  fiber/matrix  interface  can  be  written  as  follows: 


From  equation  (33): 


_  _  *o  (  Ef ) 1/2 


For  graphite/epoxy: 


0.027 

f  241  )1/a 

0.1524 

2.5  K, 

11.21/ 

K. 

where  K,  can  be  obtained  from  equation  (B-1)  and  Tables  B-1  and  B-2. 


(B-2) 


(B*3) 


(B-4) 


(B-5) 


Table  B-3.  G, 

nf/Qnf  vs- 

Gm/G,  for  Various  t,/d,  (tm/df=0.l). 

t,/d. 

=  0.001 

=0.005 

t./d, 

=  0.01 

Gr  :'G, 

& 

i-g+g2-g3 

GrVG, 

2 

i-g+gz-g3 

G^ 

2 

1-g  +  g2-g: 

1 

0.092 

0.9157 

1 

0.092 

0.9157 

1 

0.092 

0.9157 

10 

0.0965 

0.912 

10 

0.112 

0.90 

10 

0.1303 

0.885 

102 

0.1325 

0.883 

10s 

0.235 

0.810 

10z 

0.3175 

0.751 

103 

0.3175 

0.751 

103 

0.693 

0.455 

103 

0.9526 

0.09 

104 

0.9526 

0.09 

2130 

1.0 

0. 

1078 

1.0 

0. 

10680 

1.0 

0. 
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Table  B-4.  GINF/CINF  vs.  Gm/G,  for  Various  yd,  (tm/d,=  1.0). 


Gm/G, 

Vd, 

_g 

=0.0C1 

i-g+g2-g3 

Gm/Gi 

|<o 

=0.005 

i-g+g2-g3 

Gm/G, 

yd,. 

g_ 

=  0.01 
i-g+g2- 

1 

0.2258 

0.8137 

1 

0.2258 

0.8137 

1 

0.2258 

0.8137 

10 

0.2278 

0.8123 

10 

0.2348 

0.8074 

10 

0.2435 

0.8014 

102 

0.2454 

0.800 

102 

0.311 

0.756 

102 

0.3772 

0.7114 

103 

0.3791 

0.71 

103 

0.72 

0.425 

103 

0.98 

0. 

104 

1.0 

0. 

2043 

1.0 

0. 

1032 

1.0 

0. 
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